Fast and High Accuracy Multigrid Solution of the Three Dimensional Poisson Equation
نویسندگان
چکیده
We employ a fourth-order compact finite difference scheme (FOS) with the multigrid algorithm to solve the three dimensional Poisson equation. We test the influence of different orderings of the grid space and different grid-transfer operators on the convergence and efficiency of our high accuracy algorithm. Fourier smoothing analysis is conducted to show that FOS has a smaller smoothing factor than the traditional second-order central difference scheme (CDS). A new method of Fourier smoothing analysis is proposed for the partially decoupled red-black Gauss–Seidel relaxation with FOS. Numerical results are given to compare the computed accuracy and the computational efficiency of FOS with multigrid against CDS with multigrid. c © 1998 Academic Press
منابع مشابه
Optimal injection operator and high order schemes for multigrid solution of 3D poisson equation
We present a multigrid solution of the three dimensional Poisson equation with a fourth order 19-point compact nite diierence scheme. Using a red-black ordering of the grid points and some geometric considerations, we derive an optimal scaled injection operator for the multigrid algorithm. Numerical computations show that this operator yields not only the smallest overall CPU time, but also the...
متن کاملHigh Accuracy Multigrid Solution of the 3d Convection-diiusion Equation
We present an explicit fourth-order compact nite diierence scheme for approximating the three dimensional convection-diiusion equation with variable coeecients. This 19-point formula is deened on a uniform cubic grid. Fourier smoothing analysis is performed to show that the smoothing factor of some relaxation techniques with our scheme is smaller than 1. We further design a parallelization-orie...
متن کاملNumerical Solution of Multidimensional Exponential Levy Equation by Block Pulse Function
The multidimensional exponential Levy equations are used to describe many stochastic phenomena such as market fluctuations. Unfortunately in practice an exact solution does not exist for these equations. This motivates us to propose a numerical solution for n-dimensional exponential Levy equations by block pulse functions. We compute the jump integral of each block pulse function and present a ...
متن کاملA cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces
We present a finite volume method for the solution of the two-dimensional Poisson equation ∇· (β(x)∇u(x)) = f(x) with variable, discontinuous coefficients and solution discontinuities on irregular domains. The method uses bilinear ansatz functions on Cartesian grids for the solution u(x) resulting in a compact nine-point stencil. The resulting linear problem has been solved with a standard mult...
متن کاملHigh accuracy multigrid solution of the 3D convection±diusion equation
We present an explicit fourth-order compact ®nite dierence scheme for approximating the three-dimensional (3D) convection±diusion equation with variable coecients. This 19-point formula is de®ned on a uniform cubic grid. Fourier smoothing analysis is performed to show that the smoothing factor of certain relaxation techniques used with the scheme is smaller than 1. We design a parallelizatio...
متن کامل